Source of Acquisition 
NASA Washington, D C. 


mu ilium i 


uni inn inn inn iiiii mu min mi 1111 mi 


At ic~ 

d2) United States Patent 

Chidester et al. 


US007075457B1 


(10) Patent No.: US 7,075,457 B1 

(45) Date of Patent: Jul. 11, 2006 


(54) ENERGY INDEX FOR AIRCRAFT 
MANEUVERS 

(75) Inventors: Thomas R. Chidester, Mountain View, 
CA (US); Robert E. Lynch, San 
Carlos, CA (US); Robert E. Lawrence, 
Los Altos, CA (US); Brett G. Amidan, 
Kennewick, WA (US); Thomas A. 
Ferryman, Richland, WA (US); 

Douglas A, Drew, Peachtree City, GA 
(US); Robert J. Ainsworth, Dunwoody, 
GA (US); Gary L. Prothero, Corvallis, 
OR (US); Timothy P. Romanowski, 
Corvallis, OR (US); Laurent Bloch, 
Dallas, TX (US); William L. Craine, 
Vancouver, WA (US); Vincent J. 
Zaccardi, Eastchester, NY (US) 

(73) Assignee: The United States of America as 

represented by the Administrator of 
the National Aeronautics and Space 
Administration (NASA), Washington, 
DC (US) 

(-* ) Notice: * Subject to any disclaimer, the term of this 
patent is extended or adjusted under 35 
U.S.C. 154(b) by 23 days. ‘ 

(21) Appl. No.: 10/956,523 

(22) Filed: Sep. 22, 2004 

(51) Int.Cl. 

G08B 23/00 (2006.01) 

(52) U.S. Cl 340/963; 340/945; 340/946; 

340/969; 340/970 


(58) Field of Classification Search 340/945, 

340/946, 963, 967, 968, 969, 970; 701/3, 
701/121; 244/181,191 

See application file for complete search history. 

(56) References Cited 

U.S. PATENT DOCUMENTS 


4,667,196 A * 5/1987 Kaul 340/954 

5,053,767 A * 10/1991 Zweifel et al 340/968 

5,079,711 A * 1/1992 Lambregts et al. 701/3 


* cited by examiner 

Primary Examiner — Van T. Trieu 

(74) Attorney, Agent, or Finn — John F. Sehipper; Robert M. 
Padilla 

(57) ABSTRACT 

Method and system for analyzing, separately or in combi- 
nation, kinetic energy and potential energy and/or their time 
derivatives, measured or estimated or computed, for an 
aircraft in approach phase or in takeoff phase, to determine 
if the aircraft is or will be put in an anomalous configuration 
in order to join a stable approach path or takeoff path. A 
reference value of kinetic energy and/or potential energy (or 
time derivatives thereof) is provided, and a comparison 
index for the estimated energy and reference energy is 
computed and compared with a normal range of index 
values for a corresponding aircraft maneuver. If the com- 
puted energy index lies outside the normal index range, this 
phase of the aircraft is identified as anomalous, non-normal 
or potentially unstable. 

19 Claims, 3 Drawing Sheets 
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ENERGY INDEX FOR AIRCRAFT 
MANEUVERS 

ORIGIN OF THE INVENTION 

5 

This invention was made, in part, by one or more employ- 
ees of the U.S. government. The U.S. government has the 
right to make, use and/or sell the invention described herein 
without pay ment of compensation therefor, including but not 
limited to payment of royalties. to 

FIELD OF THE INVENTION 

This invention relates to monitoring and analysis of 
kinetic energy and potential energy of an ascending or 15 
descending aircraft. 

BACKGROUND OF THE INVENTION 

An aircraft that is ascending following takeoff or descend- ' 20 
ing on approach will have measurable kinetic energy and 
potential energy components, and these components will 
change with time in measurable manners. Desirable energy 
states for both takeoff and landing can be determined from 
aircraft manufacturer guidance for these phases of flight. For 25 
example, where the approach occurs at an airport with an 
operable and reliable instrument landing system (ILS), the 
ILS system may provide- data recorded on the aircraft to 
serve as a standard for comparing observed kinetic and 
potential energy components for an aircraft near the ground, 30 
below 2500 feet altitude and for an assumed straight path to 
a touchdown site. If the airport has no operable and reliable 
ILS, or if the aircraft is not near the ground, another 
mechanism for providing a standard for measurements or 
estimates is needed. On takeoff, where no electronic guid- 35 
ance comparable to the glideslope is available, the aircraft 
climb profile can be compared to manufacturer guidance or 
to observed performance for recorded aircraft departures 
from the particular airport. 

The airline industry has become concerned with the 40 
problem of unstable aircraft approaches, because approach 
and landing accidents often begin as unstable approaches. 

An “unstable approach” is often defined as an approach 
where below a threshold altitude (1000 feet for IFR and 500 
feet for VFR), the aircraft is not established on a proper glide 45 
path and with a proper air speed, with a stable descent rate 
and engine power setting, and with a proper landing con- 
figuration (landing gear and flaps extended). Airlines have 
developed approach procedures that call for abandonment of 
an approach that is determined to be unstable. 50 

Development and testing of methods for detecting atypi- 
cal flights by N.A.S.A. has revealed that high energy during 
an arrival phase (below 10,000 feet but before beginning an 
approach) is the most common reason for a flight to be 
identified as atypical or out of a statistically normal range. 55 
An atypical high energy arrival phase often corresponds to 
aircraft kinetic energy and/or potential energy that requires 
dissipation of 10-30 percent more energy than is required 
for a normal arrival phase. A normal arrival phase may 
correspond to about a 3° slope glide path and decelerating to 60 
an airspeed of about 250 knots during descent through 
10,000 feet altitude to a standard reference speed around 
2,500 feet altitude, when beginning an approach. 

More than half of the high energy arrivals identified by 
atypicality analysis were brought under control within sta- 65 
bilized approach criteria; some of the remainder of the high 
energy arrivals were abandoned. In contrast, where these 
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findings were used to define and search for a high-energy 
arrival exceedance, about three times as many excedances 
were detected; and the resulting unstable approaches were 
found to occur more frequendy than the recoveries. 

It may be possible to identify a first class of high energy 
arrivals where recovery and subsequent stabilization is pos- 
sible and relatively easy, and a second class of high energy 
arrivals in which recovery and subsequent stabilization is 
likely to be difficult or impossible. However, the present 
procedures for determining presence of a normal approach 
include an electronic glide slope that extends linearly from 
the end of a target runway to the aircraft, whereas a typical 
(normal) aircraft approach path is curved and follows the 
electronic glide slope only from about 1,800 feet above the 
field to the end of the runway. 

A 3-to-l glide path slope, corresponding to decrease of 
1 ,000 feet in altitude for every 3 nautical miles horizontal * 
travel, is often desirable during an arrival phase. The air 
speed is 250 knots or less by regulation below 10,000 feet, 
and the aircraft decelerates to a reference speed before 
joining the approach path. These parameters are directly 
available but are unlikely to prove to be the only relevant 
parameters in determining whether a flight arrival phase is 
normal or other than normal. 

What is needed is a system that (1) provides estimates of 
kinetic energy and potential energy components, and rates of 
change thereof, of an ascending or descending aircraft at any 
altitude, (2) provides reference values of kinetic energy and 
potential energy components, and rates of change thereof, of 
the aircraft; (3) provides one or more comparison indices for 
the estimated and reference values; and (4) advises an 
aircraft operator if the measured comparison indices are too 
far from the corresponding index values for a normal flight 
approach. * 

SUMMARY OF THE INVENTION 

These needs are met by the invention, which provides a 
method and a system for monitoring an ascending or 
descending aircraft to determine if the kinetic and/or poten- 
tial energy of the aircraft is within, or is outside of, a range 
for a normal flight. This invention can be used in post-flight 
review of flight data and/or as part of a flight operations 
quality assurance program, and can be implemented in an 
aircraft flight management computer to alert a pilot in real 
time to presence of an anomalous energy state. 

In one embodiment, the method includes the steps of: 

(i) providing an estimate or measurement of a value 
(referred to as an “estimated value”) of an energy compo- 
nent, E(t„)=dl*KE(t, I )+d2‘PE(t II ) of an aircraft during an 
ascent phase or descent phase of a flight, at each of a first 
sequence of times (n~l . . . , N; N^ 2), where d l and d2 are 
selected real numbers, not both 0, such as (dl,d2)=(l,0) or 
(1,1) or (0,1) or (d,l-d) with 0<d<l. 

(ii) providing a reference value of the energy component 
E(t' n ;ref) and/or reference values of the separate kinetic and 
potential energy components, KE(t' ;l ;ref) and PE(t’ ;i ;ref) at a 
time, t=t’ m determined with reference to the time t„ 
(n=l, . . . N); 

(iii) computing an index of comparison value Cl{E(t„), 
E(t 1 „;ref)} of the estimated and reference energy components 
for at least one time numbered n; and 

(iv) when the comparison index value Cl lies outside a 
selected range, interpreting this condition as indicating that 
the estimated energy component is anomalous. 

A comparison index may be based on one or more point 
values of the estimated -and reference values of the energy 
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component; or may be based on a weighted average over 
time of the estimated and reference values of the energy 
components. The comparison index may be chosen from a 
group of such indices that includes: (1) a first ratio E(t„)/E 
(t' n ;ref)}; (2) a second ratio E<t*„;ref)/E(t n ); (3) a difference 5 
E(t„)-E(t’, z ;ref)}; (4) an absolute difference (E(t„)-E(t'„;ref)l; 

(5) a normalized difference {E(t w )-E(t , /I ;ref)}/{a*E(t rt )+(l- 
a)-E(t , /J ;ref)}, where a is a selected value in a range O^a^l; 
and (6) a weighted average of the differences KE(t„)-KE 
(t’„;ref) and/or PE(t„)“PE(t , n ;ref). An analogous procedure 10 
can be applied to time rates of change of the kinetic and 
potential energy components. 

One energy index is a ratio of actual aircraft energy 
divided by ideal aircraft total energy during an arrival phase. 

If this ratio lies near a boundary but outside a “normal” 15 
range for the energy index (e.g., between about 0.90 and 
1.10), this arrival phase may be considered non-normal, and 
appropriate remedial actions may be taken to recover to a 
stabilized approach. If this ratio is below a first threshold. ‘ 
(e.g., below 0.85) or above a second threshold (e.g., above 20 
1.20) for an arrival phase, the aircraft is unlikely to be able 
to recover, and the aircraft is better advised to abandon the 
approach, to execute a go-around, and to re-enter a new 
arrival phase. 

25 

BRIEF DESCRIPTION OF THE DRAWINGS 


FIGS. 1A and IB illustrate environments in which the 
invention can be practiced. 

FIGS. 2 and 3 are flow charts of procedures for practicing 30 
an embodiment of the invention. 


DESCRIPTION OF BEST MODES OF THE 
INVENTION 


FIGS. 1A and IB illustrate environments for an ascending ■ 
aircraft (1A) and for a descending aircraft (IB) where the 
invention can be practiced. In FIG. 1A, an aircraft 11A is 
ascending, either after takeoff or in moving from a first flight 
altitude to a second flight altitude. The aircraft has an 40 
associated kinetic energy component KE(Q, measured or 
estimated or otherwise provided, at each of a first sequence 
{t n } n of two or more, time values, and has an associated 
potential energy component PE(t„), measured or estimated 
or otherwise provided, at the first sequence {t n }„ of time 45 
values. The aircraft kinetic energy and potential energy 
components are, respectively. 


KE(t)=m(t)-v( 0 2 /2+o)-/-qV2 , 


(1) 


PE(t)=m(t>g~h(t\ (2) 50 

where m(t) is the instantaneous mass (taking account of fuel 
consumption), I(t) is an instantaneous moment of inertia 
tensor for the aircraft, co(t) is an aircraft rotation vector, 
computed with reference to a center of gravity or other 55 
selected location determined with reference to the aircraft, 
(optional), v(t)=dx/dt is the instantaneous aircraft velocity 
and h(t) is the instantaneous height of aircraft eg above local 
ground. 

FIG. 2 is a flow chart of a procedure for practicing an 60 
embodiment of the invention. In step 21, an aircraft system 
measures or estimates or otherwise provides a value (re- 
ferred to as an “estimated value” for convenience herein) 
E(t n )=dl*KE(t„)+d2*PE(t n ) of an energy component of an 
aircraft during an ascent phase or descent phase of a flight, 65 
at each of a first sequence of times (n=l, . . . , Nl; N1^2), 
where dl and d2 are selected real values, not both 0. In step 
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23 , the system provides or computes a reference value 
E(t’„;ref) of the energy component at a time, t=t'„, deter- 
mined with reference to the time t n (n=l, .... . , Nl). The time 
sequence {t'„} may substantially coincide with the sequence 
{t„}, or each time value t’„ may be displaced by a calculable 
or measurable amount from the corresponding time value t„. 
In step 25 , the system computes an index of comparison 
value Cl{E(tJ, E(t’ n ;ref)} of the estimated and reference 
energy components for at least one time value pair (t„,t’„). 
When the comparison index value Cl lies outside a selected 
range for this index, the system interprets this condition as 
indicating that the estimated energy component is anoma- 
lous or non-normal or may lead to an unstable aircraft 
maneuver (step 27 ). 

A variety of comparison indices Cl can be used here. 
Some examples are: (1) a first ratio E(t n )/E(t' n ;ref); (2) a 
second ratio E(t’ n ;ref)/E(t n ); (3) a difference 
E(t„)-E(t’, z ;ref)}; (4) an absolute difference IE(t„)-E(t’„;ref)l; 
(5) a normalized difference {E(t„)-E(t , „;ref)}/{a"E(t„)+(l- 
a)*E(t’ n ;ref)}, where a is a selected real value in a range 
0~ia^l; (6) a weighted average of the differences KE(t„)- 
KE(t’ n ;ref) and/or PE(t n )-PE(t , n ;ref), such as 


w (3) 

WA = w n {\KE(tn) - me ref)\ p + \PE{t n ) - PE(t ' n ; ref)f} 1/p , 

.«= 1 


where p is a selected positive number (e.g., p=l or 2) and 
{w } n is a sequence of weight values (preferably, but not 
necessarily, non-negative); and (7) a monotonic function of 
one or more of the preceding combinations. 

The comparison index Cl may use one or a few point 
values, E(t„) and E(t'„;ref), or may use a weighted average 
of these values, such as the average 


Nl (4) 

WA1 = 2 w n Cl{E(t„), E(C, ref)). 

n= 1 


The analysis may be extended, to consider time rates of 
change, dE(t)/dt and/or dPE(t)/dt, of the energy components 
at a sequence of one or more times {t n }„ (n=l . . . , N2; 
N2^1), plus corresponding reference time rates of change, 
dE(t;ref), at a sequence of times {t’ n } n determined with 
reference to the time sequence (t n } n . Another comparison 
index, C2{dE(t„)/dt, dE(t’„;ref)/dt}, which may be the same 
or different from the comparison index Cl, is computed and 
compared with a second selected range to determine if the 
aircraft flight is anomalous or non-normal or is within a 
normal range. Again, the comparison index C2 may use 
point values or a weighted average of the values dE(t)/dt 
and/or dE(t;ref)/dt. 

FIG. 3 is a flow chart of a procedure for practicing an 
embodiment using time derivatives of the energy component 
E(t„). In step 31 , an aircraft system measures or estimates or 
otherwise provides an “estimated value”, (d/dt)E(t„)=d3*(d/ 
dt)KE(t„)+d4*(d/dt)PE(t„) of an energy component of an 
aircraft during an ascent phase or descent phase of a flight, 
at each of a first sequence of times (n=l, . . . , Nl; Nl§2), 
where d3 and d4 are selected real values, not both 0. In step 
33 , the system provides or computes a reference value 
(d/dt)E(t’ n ;ref) of the energy component at a time, t=t’„, 
determined with reference to the time t„ (n=l, Nl). The time 
sequence {t"„} may substantially coincide with the sequence 
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{t„}, or each time value t" n may be displaced by a calculable 
or measurable amount from the corresponding time value t n . 

In step 35, the system computes an index of comparison 
value C2{(d/dt)E(t„), (d/dt)E(t'„;ref)} of the estimated and 
reference energy component time deri vati ves for at least one 5 
time value pair (t„,t"„). When the comparison index value 
€2 lies outside a selected range for this index, the system 
interprets this condition as indicating that the estimated 
energy component time derivative is anomalous or non- 
normal or may lead to an unstable aircraft maneuver (step 10 
37 ). 

The analysis may be further extended to consider a third 
comparison index, C3{E(tj, E(t’ w ;ref), (d/dt)E(t„), (d/dt)E 
(t" n ;ref)}, that depends upon some or all of the estimated 
values and time rates of change of the estimated values of 15 
the energy components. Again, the comparison index C3 
may use point values or a weighted average of the values 
E(t) and/or E(t;ref) and/or (d/dt)E(t)/dt and/or (d/dt)E(t;ref). 

A formulation of, and use of, the equations of motion of 
an aircraft, including the effects of gravity, variable wind 20 
speeds, drag and lift forces on various control surfaces, 
variation of aircraft mass due to fuel consumption, and 
variable thrust, is set forth in an Appendix. A thrust vector 
is determined, as a function of the location coordinates, that 
will move the aircraft from an initial velocity vector v 0 (x 0 , 25 
y 0 ,z 0 ) to a desired final velocity vector Vj(x f , y^zj) as part of 
a takeoff phase or as part of an approach phase for a flight. 

The aircraft kinetic energy is a sum 

KE=mM 2 /2-Ko-/-cc/2, (5) .^q 

where I(t) is an instantaneous moment of inertia tensor for 
the aircraft and co(t) is an* aircraft rotation vector, computed 
with reference to a center of gravity or other selected 
location determined with reference to the aircraft. The 35 
rotational component of kinetic energy may be negligible or 
may^ be ignored for other reasons. The aircraft potential 
energy may be taken to be m-g-h, as in Eq. (2), where h is 
the aircraft altitude above local ground level. 

40 

Appendix: A Formulation of Aircraft Equations of 
Motion 

Let v(x,y,z)-(v x , v v , v,) by the vector velocity for an 
aircraft eg, let m=m(t) be the mass of the aircraft (with or 45 
without fuel consumption accounted for), and let F(x,y,z)= 

(F a ., F v , F,) by the total force vector acting upon the aircraft, 
including a thrust vector. One vector equation of motion for 
the aircraft is 

d{mv}/dt=F. (A-l) 50 

The time variation of the aircraft mass is likely close to 
linear (m(t)^mO-ml-t). When the aircraft thrust vector 
F(thrust) is known as a function of the location coordinates 
(x,y,z). Equation (A-l) can be re-expressed as 55 

d(mv)/dt = m{(8(mv) f 6x)(dx fdt) + (8[mv)/8y)(dy/d t) + (A- 2) 

(6(mv)}8z)(dzjdt)) 

oU 

= [8(mv)/8x)v x + 8(mv) f 8y)(v v + 8(mv)/8z)v z \ 

= (v / A)(mv) = F 


The force vector may be approximated as 

F(x,>’.z)=-«^^+^ r (drag/wind)+^ r (lifttK/ r (thrust) (A-3) 
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where (i,j,k) is a Cartesian coordinate system unit vector 
triad, k has the direction of a local radial vector (planar Earth 
approximation for a relatively small region), and F(drag/ 
wind) and F(lift) are a drag/wind force vector and a wind 
force vector acting on the aircraft and its control surfaces. 

A component of drag/wind force is initially assumed to 
involve only the component of total velocity in that coor- 
dinate direction so that 

F(drag/wind)=(Of J[ (v x +M J 3 2 , a y (v v +« y ) 2 , c r(v,+u 2 ) 2 ), (A-4) 

where u^u^, u v , u z ) is the local wind velocity vector, which 
may be constant or may depend upon one or more of the 
location coordinates (x,y,z). The drag coefficient vector 
o=(0 x , G y , G z ) will depend upon the angle of attack % upon 
a'ileron surfaces orientation angle a, upon flap angle/exten- 
sion p, upon rudder angle y, upon elevator angle 8, and upon 
angular orientation angles (((>,0) of the aircraft fuselage and 
empennage relative to the local Cartesian coordinate system. 
As a first approximation, the drag coefficient vector is 
expressed as a sum of terms 

o=a(aileron)+a(flap)+<7(rudder)+a(elevator)+o(fuse- 

lage)+a(enipennage), (A-5) 

where the individual vector contributions (e.g., o(aileron)) 
are determined for the particular aircraft configuration and 
projected area of the relevant control surface, using empiri- 
cal and/or experimental information. For example, if a 
particular aircraft control surface is planar and aircraft 
velocity is subsonic, a first approximation to the drag force 
coefficient for an airfoil surface can be expressed in terms of 
momentum transfer rate as 

<j(T)=alsin 2 ^, (A-6) 

where T is an angle of the airfoil surface normal relative to 
a vector representing movement of air past the airfoil. 

The lift component F(lift) depends upon total aircraft 
velocity, (v+u) 2 , upon angle of attack x of a wing or other 
contributing surface that contributes to lift, and upon local 
air density p, 

F(lift)=F«v+w) 2 , X,p) (A-l) 

preferably using Bernoulli’s equation. Unlike the drag force, 
a movement of air in one direction may give rise to a lift 
force in a different (e.g., perpendicular) direction. 

Incorporating these characterizations, the equations of 
motion of the aircraft in response to extraneous forces (wind, 
draft, lift, gravity, etc.) is re-expressed as 

( v- A)(mv)=/{ 2o x (v x +w x ) 2 +2F x (lift) } 2F x (thrust) }+y { 2<5 y 
( v v +« v ) 2 +2F v (lift) }+2F v ( thrust) }+&{-2 mg+2a 2 
( v,4-M z ) 2 +2F z (lift) }+2F 2 (thrust) } . ( A-8) 

Beginning at an initial location, (x 0 ,y 0 ,z 0 ), one goal is to 
bring the aircraft from an initial condition v(x 0 ,y 0 ,z 0 ) to. a 
desired final condition V(x / ,y / ,z / ) as the aircraft enters an 
approach zone or leaves a takeoff zone, by prescribing a 
thrust vector F(thrust) that will move the aircraft from the 
initial velocity condition v(x 0 ,y 0 ,z 0 ) to the final velocity 
condition v(x / ,yy > z / ) without violating limits on the aircraft 
variables. Other formulations of aircraft equations of motion 
can be provided that incorporate the effects of (1) variation 
of aircraft mass, (2) presence of (variable) wind velocity, (3) 
components of drag force induced on various control sur- 
faces and (4) gravity, as well as other effects that have 
smaller influence on aircraft motion. However, these equa- 
tions of motion appear to manifest the major features, 
including the effects of a programmed thrust vector to move 
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the aircraft from an initial condition v 0 (x 0 ,y 0 ,z 0 ) to a desired 
final condition vyt'x^y^Zf) as part of an ascent phase or 
descent phase of the flight. 

As a first approach, parameters ' and associated forces 
(drag, lift, thrust, etc.) may be measured at a sequence of 5 
times for a moving aircraft, and Eq. (A-8) may be solved for 
a sequence of corresponding location coordinates (x,y,z) to 
determine or estimate an aircraft velocity vector v and 
altitude h to determine aircraft kinetic energy KE and/or 
potential energy PE for this particular flight movement. 10 
Here, the kinetic energy and/or potential energy of the 
aircraft are determined or estimated off-line, after this por- 
tion of the flight is completed. 

As a second approach, a proposed aircraft maneuver to 
move from an initial velocity condition v(x 0 ,y 0 ,z 0 ) and 15 
initial altitude to a final velocity condition v(xQ,y 0 ,z 0 ) and 
final altitude can be posited and a thrust field F(thrust) can 
be determined that will accomplish this can be determined, 
through solution of Eq. (A-8), before or during execution of 
the aircraft maneuver (on-line). 20 

What is claimed is: 

1. A method of monitoring energy components of an 
aircraft in flight, the method comprising: 

providing an estimate or measurement of a value (referred 
to as an “estimated value’') of an energy component, 25 
E(t„)=dl-KE(t„)+d2-PE(t n ), of a combination of a 
kinetic energy component KE(tJ and a potential 
energy component PE(t„) of an aircraft during an ascent 
phase of a flight, at each of a first sequence of times 
(n=l, . . . , N; N^2), where dl and d2 are selected real 30 
numbers, at least one being non-zero; 
providing a reference value E(t' n ;ref) of the energy E(t n ) 
at a time, t=t’ n , determined with reference to the time t n 
(n=l, . . . , N); 

providing an index of comparison Cl {E(t n ), E(t'„;ref)} of 35 
the estimated and reference energy components for at 
least one time numbered n; and 
when a value of the comparison index Cl lies outside a 
selected range, interpreting this condition as indicating 
that the estimated energy component is anomalous. 40 

2. The method of claim 1, further comprising: 
providing reference values, KE(t’„;ref) and PE(t’ rt ;ref), of 

kinetic energy and potential energy components for 
said estimated values; and 

choosing said comparison index from the group of indices 45 
consisting of: (1) a first ratio E(t J1 )/E(t’„;ref)}; (2) a 
second ratio E(t' n ;ref)/E(t„); (3) a difference E(t„)-E 
(t’ n ;ref)}; (4) an absolute difference IE(t„)-E(t’ n ;ref)l; 

(5) a normalized difference {E(t n )-E(t , w ;ref)}/{a-E(t„)+ 
(l-a)E(tV,ref)}, where a is a selected value in a range 50 
0^a^ 1 ; and (6) weighted averages, 
IKE(t n )-KE(t’„;ref)i p and iPE(t n )-PE(t , „;ref)l p , of differ- 
ences of kinetic energy terms and of potential energy 
terms, where p is a selected positive number. 

3. The method of claim 1, further comprising choosing 55 
said comparison index to include at least one of said values 
E(t„) for said estimated energy component and at least one 

of said values E(t n .;ref) for said reference energy component. 

4. The method of claim 1, further comprising choosing 
said comparison index to include at least one weighted 60 
average of said values E(t„) for said estimated energy 
component and at least one weighted average of said values . 
E(t„.;ref) for said reference energy component, over said 
respective sequences of times {t„} and {t’ n }. 

5. The method of claim 1, further comprising: 65 

providing an estimated value (d/dt)E(t n )=d3*(d/dt)KE 

(t„)+d4-(d/dt)PE(t n ) of a time rate of change of said 
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estimated energy component, where d3 and d4 are 
selected real values, not both 0; and 
providing a reference value (d/dt)E(t" n ;ref) of a time rate 
of change of said reference energy component, where 
t" n is a time determined with reference to said time t„. 

6. The method of claim 5, further comprising: 
providing an index of comparison C2{(d/dt)E(t„); (d/dt) 

E(t M „ ;ref)} of said time rates of change of said esti- 
mated energy component and said reference energy 
component for at least one time numbered a”. 

7. The method of claim 6, further comprising choosing 
said comparison index to include at least one of said values 
(d/dt)E(t n ) for said estimated energy component and at least 
one of said values (d/dt)E(t" n ;ref) for said reference energy 
component. 

8. The method of claim 6, further comprising choosing 
said comparison index to include at least one weighted 
average of said values (d/dt)E(t w ) for said estimated energy 
component and at least one weighted average of said values 
(d/dt)E(t" n ;ref) for said reference energy component, over 
said respective sequences of times {t n -} and {t" n .}. 

9. The method of claim 6, further comprising: 

when a value of said comparison index C2 lies outside a 
selected range, interpreting this condition as indicating 
that said time rate of change of said estimated energy 
component is anomalous. 

10. The method of claim 5, further comprising: 
providing an index of comparison C3{E(t n .), E(t w .;ref); 

(d/dt)E(t„), (d/dt)E(t"„;ref)} of said estimated values 
and time rates of change of said estimated energy 
component and said reference energy component for at 
least one time numbered n'; and 
when a value of the comparison index C3 lies outside a 
selected range, interpreting this condition as indicating 
that at least one of said value and said time rate of 
change of said estimated kinetic energy component is 
anomalous. 

11. The method of claim 1, further comprising choosing 
said real numbers dl and d2 to be one of the following pairs 
of real numbers: (dl,d2)=(l,0), (dl,d2)=(l,l), (dl,d2)=(0,l), 
and (dl,d2)=(d,l-d) with 0<d<l. 

12. The method of claim 1, further comprising determin- 
ing said energy component by a process comprising: 

measuring aircraft thrust vector components and at least 
one of extraneous force vector components, comprising 
wind vector components, drag force vector compo- 
nents, lift force vector components and gravity force 
components, at said sequence of times {t„}; 
determining or estimating at least one of velocity vector 
v of said aircraft and altitude of said aircraft at location 
coordinates (x n ,y n ,z 7I ) corresponding to said time t n , 
from a solution of an equation (v-A)(mv)=a sum of 
extraneous vector forces F(extraneous)acting on said 
aircraft+a thrust vector force F(thrust).provided by one 
or more engines of said aircraft; 
determining or estimating at least one of an aircraft 
velocity vector v(x rt ,y„,z n ) and an aircraft altitude for at 
least one of said times t„; and 
computing said estimated or measured value of said 
energy component E(t„) for at least one of said times t n , 
using the determined or estimated aircraft velocity 
vector and aircraft altitude. 

13. The method of claim 1, further comprising determin- 
ing said energy component by a process comprising: 

providing at least one of extraneous force vector compo- 
nents F(extraneous), comprising wind vector compo- 
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nents, drag force vector components, lift force vector 
components and gravity force components, at said 
sequence of times {t„}; 

providing an estimate of thrust vector components 
F(thrust) required to transport an aircraft from a 
selected initial velocity condition v(x 0 ,y 0 ,z 0 ) to a 
selected final velocity condition vfx^y. y,z r ) under influ- 
ence of the at least one extraneous force components; 
determining or estimating at least one of velocity vector 
v of said aircraft and altitude of said aircraft at location 
coordinates (x„,y„,z n ) corresponding to said time t w , 
from a solution of an equation (v-A)(mv)=a sum of 
extraneous vector forces F(extraneous)acting on said 
aircraft+the thrust vector force F(thrust), at said 
sequence of times {t rt }; 

determining or estimating at least one of an aircraft 
velocity vector v(x n ,y„,z„) and an aircraft altitude for at 
least one of said times t w ; and 
computing said estimated or measured value of said 
energy component E(t n ) for at least one of said times t„, 
using the determined or estimated aircraft velocity 
vector and aircraft altitude. 

14. A method of monitoring energy components of an 
aircraft in flight, the method comprising: 

providing an estimate or measurement of a value (referred 
to as an “estimated value”) of an energy component, 
(d/dt)E(t„)==d3*(d/dt)KE(t M )4d4-(d/dt)PE(t„), of a com- 
bination of time derivatives of a kinetic energy com- 
ponent KE(t rt ) and a potential energy component PE(t n ) 
of an aircraft during an ascent phase of a flight, at each 
of a first sequence of times (n=l, . . . , N; N^2), where 
d3 and d4 are selected real numbers, at least one being 
non-zero; 

providing a reference value (d/dt)E(t'* w ;ref) of the energy 
component time derivative (d/dt)E(t n ) at a time, t=t"„, 

determined with reference to the time t ?J (n=l N); 

providing an index of comparison C2{(d/dt)E(t„), (d/dt) 
E(t”„;ref)} of the estimated and reference energy com- 
ponent time derivatives for at least one time numbered 
n; and 

when a value of the comparison index C2 lies outside a 
selected range, interpreting this condition as indicating 
that the estimated energy component is anomalous. * 

15. The method of claim 14, further comprising: 
providing reference values, (d/dt)KE(t’ n ;ref) and (d/dt)PE 

(t'„;ref), of kinetic energy and potential energy compo- 
nent time derivatives for said estimated values; and 
choosing said comparison index from the group of indices 
consisting of: (1) a first ratio 

(d/dt)E(t n )/(d/dt)E(t M „;ref)}; (2) a second ratio d/dt)E 
(t H rt ;ref)/(d/dt)E(t n ); (3) a difference d/dt)E(tJ-(d/dt)E 
(t’Vjef); (4) an absolute difference ld/dt)E(t n )"(d/dt)E 
(t" n ;ref)l; (5) a normalized difference {d/dt)E(t„)-(d/dt) 
E(t" n ;ref) }/ { a- .d/ dt)E(t n )+( 1 -aHd/ dt)E(t"„ ;ref) } , 
where a is a selected value in a range O^a^ 1; and (6) 
weighted averages, l(d/dt)KE(t ?I )-(d/dt)KE(t" w ;ref)! p 
and !(d/dt)PE(t M )-(d/dt)PE(t" n ;ref)! p , of differences of 
kinetic energy terms and of potential energy terms, 
where p is a selected positive number. 
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16. The method of claim 14, further comprising choosing 
said comparison index to include at least one weighted 
average of said values (d/dt)E(t„) for said estimated energy 
component and at least one weighted average of said values 

5 (d/dt)E(tV;ref) for said reference energy component, over 
said respective sequences of times {t,J and {t'„}. 

17. The method of claim 14, further comprising choosing 
said real numbers dl and d2 to be one of the following pairs 
of real numbers: (d3,d4)=(l,0), (d3,d4)=(l,l), (d3,d4)=(0,l), 

10 and (d3,d4)=(d,l-d) with 0<d<l. 

18. The method of claim 14, further comprising deter- 
mining said energy component by a process comprising: 

measuring aircraft thrust vector components and at least 
one of extraneous force vector components, comprising 
15 wind vector components, drag force vector compo- 
nents, lift force vector components and gravity force 
components, at said sequence of times {t n }; 

determining or estimating at least one of velocity vector 
v of said aircraft and altitude of said aircraft at location 
20 coordinates (x„,y„,z n ) corresponding to said time t„, 
from a solution of an equation (v*A)(mv)=a sum of 
extraneous vector forces F(extraneous)acting on said 
aircraft+a thrust vector force F(thrust) provided by one 
or more engines of said aircraft; 

25 determining or estimating at least one of an aircraft 
velocity vector v(x„,y„,z n ) and an aircraft altitude for at 
least one of said times t n ; and 

computing said estimated or measured value of said 
energy component (d/dt)E(t n ) for at least one of said 
30 times t„, using the determined or estimated aircraft 
velocity vector and aircraft altitude. 

19. .The method of claim 14, further comprising deter- 
mining said energy component by a process comprising: 

providing at least orie of extraneous force vector compo- 
35 nents F(extraneous), comprising wind vector compo- 
nents, drag force vector components, lift force vector 
components and gravity force components, at said 
sequence of times {t„}; 

providing an estimate of thrust vector components 
40 F(thrust) required to transport an aircraft from a 
selected initial velocity condition v(x 0 ,y 0 ,z 0 ) to a 
selected final velocity condition VCx^y^z^) under influ- 
ence of the at least one extraneous force components; 

‘ determining or estimating at least one of velocity vector 
45 v of said aircraft and altitude of said aircraft at location 
coordinates (x„,y n ,z n ) corresponding to said time t n , 
from a solution of an equation (v-A)(mv)=a sum of 
extraneous vector forces F(extraneous)acting on said 
aircraft+the thrust vector force F(thrust), at said 
50 sequence of times {t„}; 

determining or estimating at least one of an aircraft 
velocity vector v(x n ,y n ,z w ) and an aircraft altitude for at 
least one of said times t„; and 

computing said estimated or measured value of said 
55 energy component (d/dt)E(t w ) for at least one of said 
times t„, using the determined or estimated aircraft 
velocity vector and aircraft altitude. 

* * $ * * 



